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ABSTRACT 

Context. The initial-final mass relationship (IFMR) for stars is important in many astrophysical fields, such as the evolution of 
galaxies, the properties of type la supernovae (SNe la) and the components of dark matter in the Galaxy. 
Aims. The purpose of this paper is to obtain the dependence of the IFMR on metallicity. 

Methods. Following Paczyhski & Ziolkowski (1968) and Han et al. (1994), we assume that the envelope of an asymptotic giant 
branch (AGB) or a first giant branch (FGB) star is lost when the binding energy of the envelope is equal to zero (AW = 0) and the 
core mass of the AGB star or the FGB star at the point (AW = 0) is taken as the final mass. Using this assumption, we calculate the 
IFMRs for stars of different metallicities. 

Results. We find that the IFMR depends strongly on the metallicity, i.e. Z = 
0.0001,0.0003,0.001,0.004,0.01,0.02,0.03,0.04,0.05,0.06,0.08 and 0.1. From Z = 0.04, the final mass of the stars with a 
given initial mass increases with increasing or decreasing metallicity. The difference of the final mass due to the metallicity may 
be up to 0.4 M G . A linear fit of the initial-final mass relationship in NGC 2099 (M37) shows a potential evidence of the effect of 
metallicity on the IFMR. The IFMR for stars of Z = 0.02 obtained in the paper matches well with those inferred observationally 
in the Galaxy. For Z > 0.02, helium WDs are obtained from the stars of M\ < 1.0M o and this result is upheld by the discovery of 
numerous low-mass WDs in NGC 6791 which is a metal-rich old open cluster. Using the IFMR for stars of Z = 0.02 obtained in the 
paper, we have reproduced the mass distribution of DA WDs in Sloan DR4 except for some ultra-massive white dwarfs. 
Conclusions. The trend that the mean mass of WDs decreases with effective temperature may originate from the increase of the initial 
metallicities of stars. We briefly discuss the potential effects of the IFMR on SNe la and at the same time, predict that metal-rich 
low-mass stars may become under-massive white dwarfs. 
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1. Introduction 

White dwarfs (WDs) are the endpoint of the evolution of stars 
with initial masses ranging from about 0.1 M Q to about 8 M Q . 
The vast majority of stars in the Galaxy belong to the mass 
range and over 97% of the stars in the Galaxy will eventually 
end up as WDs (Fontaine et al. 2001 ). As a result, WDs may 
give direct information about star formation during the Galaxy's 
earliest epochs (Gates et al. 2004). WDs show their importance 
in many fields. It has been about 50 yrs since Schmidt (1959) 
recognized the usefulness of WDs as cosmochronometers. WDs 
may also be the dominating component of dark matter in the 
Galaxy (Alcock et al. 1 19991 Chabrier [19991 . As a matter of fact, 
WDs are very important for type la supernovae (SNe la) since 
it is believed that SNe la are from the thermonuclear runaway 
of carbon-oxygen white dwarfs (CO WD) (see the reviews by 
Hillebrandt & Niemeyer (|2000T > and Leibundgut (l2000l i). It is 
well known that the masses of WDs are typically of the or- 
der of half that of the Sun, while their radii are similar to that 
of a planet. However, a detailed knowledge of WDs is still un- 
clear in observation and theory (Fontaine et al. 1200 11 Moroni & 
Straniero 2002, 2007). Among all the uncertainties of WDs, the 
correlation between initial-final mass relationship (IFMR) and 
metallicitiy is a very important one. It is well known that low 
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metallicity leads to a larger CO WD for a given initial mass with 
Z < 0.02 (Umeda et al. 1999a). However, it is necessary to check 
the cases of Z > 0.02 since Z € [0.06 - 0.1] is possible for some 
ultra-luminous galaxies (Roberts & Hynes 1994; Ruiz-Lapuente 
et al. [19951 Terlevich & Forbes [20021). 

The IFMR for stars over a large mass range (e.g. 0.8-8 M Q ) 
is a powerful input to chemical evolution models of galaxies 
(including enrichment in the interstellar medium) and there- 
fore can enhance our understanding about star formation effi- 
ciencies in these systems (Ferrario et al. 2005, Kalirai et al. 
2007b ). The IFMR is also an important input for modelling the 
luminosity functions of Galactic disk WDs and the cooling se- 
quences of halo clusters, which may directly yield the age of 
the Galactic disk and halo components (Kalirai et al. 2007b). 
Meanwhile, the IFMR represents the mass loss of a star over its 
entire evolution and it is possible to get some indications about 
the origin and evolution of hot gas in elliptical galaxies from the 
IFMR (Mathews 119901 ). Since Weidemann d 19771 ) showed the 
first comparison between observations and theoretical predica- 
tions of the IFMR, many observations gave constraints on the 
relationship by studying the properties of white dwarfs in open 
clusters or field white dwarfs, especially during the last decade 
(Herwi g [19951 Reid [19961 Koest er & Reimers [T9961 F inley & 
Koes ter [19971 Claver et al |200T1 Williams et al. 120041 Ferrario 
et al. 120051 Williams et al. 120061 ). Some empirical relations have 
also been suggested based on different observations (Weidemann 
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[19831120001 Ferrario et al. 120051 Williams et al.[2006 Dobbie et 
al. [2006b). At the same time, numerous sets of stellar evolution 
models have been calculated to study the relationship in theory 
and to give the IFMRs for stars of different metallicity (Han et 
al. 1 19941 Girardi et al. 2000). Although great progress has been 
made in observation and theory, there still exist many uncertain- 
ties, i.e. what is the origin of the intrinsic scatter of WD mass or 
whether there is any dependence of the IFMR on the metallic- 
ity (Kalirai et al. 2005; Williams 2006 ). The correlation between 
the IFMR and the metallicity has been established in theory for 
many years (Han et al. 119941 Girardi et al. 120001) . but the evi- 
dence of the dependence of the IFMR on the metallicity was not 
found until 2005 (Kalirai et al. 120051 ) . Many theoretical calcula- 
tions showed that the core mass at the first thermal pulse (TP) 
in the asymptotic giant branch (AGB) may be taken as the final 
mass (see the review by Weidemann 2000). However, these the- 
oretical studies only focused on some special metallicities and 
it is also difficult for some stellar evolution codes to determine 
which is the first thermal pulse. In the paper, we use a simple 
but robust method to systemically study the IFMR over a wide 
metallicity range, i.e. 0.0001 < Z < 0.1. We also briefly discuss 
the potential effect of the IFMR on SNe la. 

In section[2] we describe our model and physical inputs. We 
give the results in section [3] and show discussions and conclu- 
sions in section |4] 

2. Model and Physical Inputs 

2.1. the model 

When a star is in the AGB stage, its envelope may be blown off if 
the binding energy (BE) of the envelope changes from negative 
to positive (Paczyhski & Ziolkowski 1968 ). It is clearly shown 
in Fig. 1 of Han et al. ( 119941 ) that the BE of an AGB star changes 
from negative to positive through the AGB evolution (also Fig. 
[2]in this paper). The BE of the envelope can be calculated by 

r M * Gm 

AW = I ( + U)dm, (1) 

Jm c r 

where M c is the core mass, M s is the surface value of the mass 
coordinate m, and U is the internal energy of thermodynamics 
(including terms due to ionization of H and dissociation of H2, 
as well as the basic j%T Ifi for a perfect gas). We assume that 
the envelope of a star is lost if AW = and the core mass at 
this point is the final WD mass. The virtue of this method is 
that the mechanism of mass loss need not be considered. Han 
et al. d 19941 ) used this assumption to get the IFMR for stars of 
Z = 0.02 and Z = 0.001. The binary population synthesis (BPS) 
incorporating their IFMRs well reproduced the mass distribution 
of planetary nebula nuclei (PNN). Their result indicated that this 
method is robust to calculate the final mass for a star of given 
initial mass and metallicity (see Han et al. 1994 for details). 

In principle, the envelope has enough energy to escape to 
infinity when AW > 0. As mentioned in Han et al. (1994), how- 
ever, it is not clear whether the instant envelope ejection occurs 
or not. Radiation might take away some of the energy which is 
required to go into outward motion in order to achieve the enve- 
lope ejection, but it is possible that a series of oscillations on a 
dynamical time-scale with amplitude growing occur until AW is 
large enough to energize envelope ejection (Han et al. 1994 ). The 
discussion above means that AW > may only be a necessary 
condition to eject the envelope of an AGB star. Our assumption 
implies that a superwind starts at or after AW = 0. Fortunately, 



the change of the core mass during the superwind phase is neg- 
ligible at low masses and relatively modest at high masses (Han 
etal. [T994l ). 

For saving cpu time, we simply treated the average evolution 
of thermally pulsing AGB models, i.e. skipping thermal pulses 
by taking a longer time-step. We give some discussions on the 
influence of TPs and mass loss in subsection |4.1| 

2.2. Physical Inputs 

We use the stellar evolution code of Eggleton (119711 119721 
119731 ), which has been updated with the latest input physics over 
the last three decades (Han et al. [19941 Pols et al. [19931 [19981 . 
We set the ratio of mixing length to local pressure scale height, 
a = l/H p , to 2.0, and set the convective overshooting parame- 
ter, 6 ov , to 0.12 (Pols et al. [19971 Schroder et al. [19971 which 
roughly corresponds to an overshooting length of 0.25//p. The 
range of metallicity is from 0.0001 to 0.1, i.e. 0.0001, 0.0003, 
0.001, 0.004, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.08, 0.1. The 
opacity tables for these metallicties are compiled by Chen & 
Tout (2007 ) from Iglesias & Rogers (1996) and Alexander & 
Ferguson ( 1994). For a given Z, the initial hydrogen mass frac- 
tion is assumed by 

X = 0.76-3.0Z, (2) 

(Pols et al. 119981 ), and then the helium mass fraction is Y = 
1 - X - Z. Based on the correlation among X, Y and Z used here, 
Pols et al. ( 1998 ) accurately reproduced the color-magnitude di- 
agrams (CMD) of some clusters. 

The mass loss might affect the IFMR. Han et al. (1994) have 
shown that the steady stellar wind based on the observational 
relation given by Judge & Stencel d 1 99 11 ) can not significantly 
affect the IFMR. We also test the effect of Reimers' wind (1975) 
with 77 = 1 / 4 on th e IFMR of Z = 0.02. A si milar re sult to that 
in Han et al. (1994) is obtained. Buzzoni et al. (2006 ) also noted 
that Reimers' wind parametrization poorly reproduce the relative 
number of planetary nebulae (PNe) in late-type galaxies, while 
a better fit for the relative number of PNe is obtained using the 
empirical IFMR in Weidemann (2000). This might imply that 
the Reimers' wind is a unreasonable mechanism to calculate the 
IFMR. We therefore do not include any wind mass loss in our 
calculation for any of the metallicities. 

3. Results 

3.1. the initial-final mass relationship for stars of different 
metallicities 

Fig. [TJ shows the IFMRs for stars of different metallicities. The 
IFMRs in the paper may be approximated by two parabola seg- 
ments, which is similar to that of Han et al. ( 1994]) for Z — 0.02 
and Z = 0.001. The fitted formulae of the IFMRs for stars of 
different metallicities are shown in Appendix A. 

Han et al. (1994 ) showed that the BE of the envelope of Pop 
I low-mass stars (M; < 1 .0) may get the point of AW = on the 
first giant branch (FGB) and helium WDs may be obtained. This 
phenomenon was used to explain the existence of low mass WD 
in Stein 205 IB and 40 EriB. We get a similar result to that of 
Han et al. (119941 ) for those stars with Z > 0.02 and Afj < l.OM 
(See those Mf < 0.5M o in Fig. [TJ. The final masses of the He 
WDs for a given initial mass slightly depend on the metallicity. 
For a given initial mass, a larger He WD is obtained for a small 
metallicity (see in Table [TJ. This is because for a perfect gas, 
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Fig. 1. Initial-final mass relationship (IFMR) for stars of different metallicities. The points of Mf < O.5M represent helium WDs. 
The thick solid lines are the mass threshold for the second dredge-up (see the text for details) and the thick dotted ones show the 
transition between carbon-oxygen white dwarfs and oxygen-neon-magnesium ones in our models. The dashed line in the lower 
panel shows the IFMR for stars of Z = 0.04 as a comparison. 



mean molecular weight, jj, increases with metallicity according 
to equation (2) (Kippenhahn & Weigert fl 9901 1 and the thermal 
energy therefore decreases with the metallicity. Though the fi- 
nal radii of the low-mass stars when AW = in FGB slightly 
increases with the metallicity (see in Table [2]), it does not sig- 
nificantly affect the gravitational energy of the envelope. So, the 
core mass dominating the gravitational energy of the envelope 
decreases with the metallicity to counteract the decrease of the 
thermal energy. 

According to the above finding, there are probably several 
under-massive WDs (Mwd < O.5M ) in metal-rich old clus- 
ters. Observationally, Kalirai et al. d2007al l recently found that 
most WDs are under-massive and the mean mass of the WDs is 
< M >= 0.43 ±O.O6M in NGC 6791, which is one of the oldest 
(8-12Gyr) and most metal-rich ([Fe/H] =+0.4+0. 1) open clusters 
in the Galaxy (Peter son & Green [T99"8l Chaboyer et al [19991 
Stetson et al. [20031 Gratton et al. [20061 Origlia et al. [2006b . 
According to the high metallicity and the old age of NGC 6791, 
the stars with M\ ~ l.OM in NGC 6791 likely have left the 



main sequence and moved into the FGB phase (the life of stars of 
Mi = 1 .OM and Z = 0.04 - 0.06 on the main sequence is about 
7.5-9 Gyr, Chen & Tout P2"0"07l l. Kalirai et al. d2007al) suggested 
that the under-massive stars are helium WDs which have expe- 
rienced enough mass loss on the FGB to avoid the helium flash, 
which naturally resolves the age discrepancy of the cluster from 
white dwarf cooling theory (2.4 Gyr) and from main-sequence 
turnoff (8 - 12 Gyr) (Bedin et d. l20l)5l 

It is obvious from FigQ]that there is a dependence of IFMR 
on the metallicity. The IFMR gets to its lower limit at Z = 0.04, 
i.e. the final mass is smallest for Z = 0.04 for a given initial 
mass. It is well known that for stars with Z < 0.02, the effect of 
decreasing metallicity is similar to that of increasing mass, i.e. 
the evolution track of a low-mass star with a low metallicity is 
similar to that of a larger one with a higher metallicty (Umeda 
et al. |1999a) . The result of Chen & Tout d2007b reconfirmed the 
fact and suggested an opposite trend when Z > 0.04, i.e. for 
the stars with a given initial mass, the age decreases with the 
metallicity and the luminosity increases with the metallicity (see 
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Fig. 2. The change of binding energy AW with the different 
choices of core mass, M c , (equation (1)) at three evolutionary 
stages for two AGB stars, bottom panel is for M\ — 2.OOM , Z; = 
0.02 and the top panel is for M\ = 4.52M ,Zi = 0.1. The 
numbers present time sequence. The dependent variable, y = 
sinrr'(AW/10 45 erg), keeps the sign of AW/10 45 erg and is a loga- 
rithmic function of AW/10 45 erg when its absolute value is large, 
or a linear function when its absolute value is small. The dashed 
lines show the boundaries of the cores of the models defined in 
this paper, while the dot-dashed ones show the positions of he- 
lium burning shells. 



Fig. 2 and Table 1 in that paper). The phenomenon that the effect 
of the metallicity on the evolution of a star is not monotonic is 
derived from the correlation among X, Y and Z used in this paper. 
The three components all contribute to opacity which dominates 
the evolution of a star with a certain mass. As is well known, the 
opacity increases with X or Z and decreases with Y. The effect 
of increasing opacity is similar to that of decreasing mass. For 
Z < 0.04, the change of Z dominates the change of the opacity, 
while the changes of Y and X dominate the change of the opacity 
forZ>0.04. 

According to our calculation, the metallicities will result in 
a scatter of the final masses. For a certain mass, the scatter may 
be up to 0.4 M Q , which depends on the differences of metallicity 
and initial stellar mass. Williams (2006 ) compared the IFMR for 
stars in NGC 2099 (M37, Z=0.01 (Kalirai et al. 12005V ) with that 
in other two clusters ( Hyades and Praesepe, Z=0.02 (Perryman 
et al. 119981 Claver et al. 12001)) and claimed that "any metallicity 
dependence of the IFMR for M\ as 3M must be smaller than 
AM f ^ 0.05 M Q ". The difference of the final mass for Mj * 3M 
between Z = 0.02 and Z = 0.01 in this paper is about O.O4M , 
which is consistent with the result of Williams (2006 ). 

We assume that the remnant is a CO WD if carbon and oxy- 
gen have not been ignited in a star when AW = (see the thick 
dotted lines in Fig. [TV In Table[3] we show the maximum masses 
of the CO WDs, , for different metallicities. Note that the 
values of the M™ x probably have some differences from the real 
ones because of the grid density of the models in this paper. In 



our study, the are 1.2431M and 1.2132M for Z = 0.08 
and Z = 0.1, respectively. Generally, the mass of a CO WD is 
less than 1.20 M Q , otherwise carbon will be ignited in the re- 
gion with the highest temperature and the final remnant after the 
ejection of its envelope is an oxygen-neon-magnesium (ONeMg) 
WD. The overestimated values in our models for Z = 0.08 and 
Z = 0.1 are mainly from the determination of the core mass 
M c . As discussed by Han et al. (119941 ), it is difficult to deter- 
mine where the boundary of the core is. As an examination, we 
choose different core masses M c from the center to the surface 
and calculated the BE AW as a function of M c . Three evolu- 
tionary stages, i.e. before, at and after AW = 0, are selected to 
check the evolution of AW = for two AGB stars. The results 
are shown in Fig. [2] In the figure, we see that there are three por- 
tions: (i) an inner region where AW increases slowly; (ii) a zone 
where AW increases sharply, the zone including the H-burning 
shell; and (iii) an outer portion where AW varies slowly with M c . 
We take the value near but outside the transition between (ii) and 
(iii) as M c (the dashed lines in the figure) to be sure that AW is 
not very sensitive to M c . This choice is good enough for stars 
with a thin transition between (ii) and (iii), i.e. both M\ and Z are 
not very large such as M\ = 2.OOM and Z = 0.02 in the figure. 
However, for those with large initial mass and high metallicity, 
the transition between (ii) and (iii) is thick (see the model of 
Mi = 4.52M ,Zi = 0.1 in Fig. [2]) and it is difficult to choose the 
boundary of the core. This may result in an uncertainty of the 
final mass by O.O3M for the model of M { = 4.52M ,Zi =0.1. 
Actually, for most of the models in the paper, our choice of M c 
is not a serious problem. The maximum uncertainty of the final 
masses derived from the choice of M c is about O.O4M , which is 
from the model of Mi = 5.50M o and Z = 0.1. 

The value of Af£° might also be overestimated if the star 
reaches AW = again after the ejection of H-rich envelope. 
There is still a helium envelope around the CO core after a 
star loses its hydrogen-rich envelope and the helium envelope 
is thicker with Mi and Z. For example, the mass of the helium 
envelope after the loss of the hydrogen-rich envelope may be as 
large as 0.2 M Q for a star of Mi = 4.52M with Z = 0.1 (see 
in Fig [2] the dot-dashed lines show the positions of the helium 
burning shells). Then, the star becomes a helium red giant. With 
the expansion of the helium envelope, the envelope might get 
the point where AW = and a part of the envelope might be lost 
again, which may lead to an overestimate of M™ x (Han et al. 
IT994V 

In Fig. [2] we see that, although the shape of AW with M c has 
not changed during the AGB phase, the time sequence is differ- 
ent for the two stars, i.e. the region of the sharp increase of AW 
moves outward with time for the low-mass star, while it is oppo- 
site for the high-mass one. This difference comes from whether 
the second dredge-up occurs or not in a star. Generally, for Pop 
I stars (Z = 0.02), the second dredge-up occurs in early-AGB 
(E AGB) stars defined by Iben (fT983l if Mi > 3.5M (Busso et 
al. 1 19991 ), which reduces the mass of H-exhausted core (see Fig. 
6 of Meng et al. 120061 ). AW reaches zero just at the process of 
core decreasing as shown in the upper panel of Fig. [2] For the 
stars with Mi < 3.5M , the second dredge-up does not occur 
and the core mass always increases with time as shown in the 
bottom panel of Fig. [2] When a star has a mass around 3.5M , 
the core mass is almost constant with the expansion of its enve- 
lope until AW = (see Fig. 6 of Meng et al. 2006). The mass 
threshold for the second dredge-up changes with metallicity Z. 
We presented the mass threshold from our models for various Z 
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in Table|4](see the thick solid line in Fig. [TJ Q From Table|4] we 
see that, since the stellar evolution is non-monotonic with Z, Q 
the mass threshold is also non-monotonic with Z. 

We notice that the final mass of the star of M\ = 5.5M with 
Z = 0.1 is larger than 1.5 M Q . It is likely that this remnant is not 
a white dwarf. When the star first arrives at the point of AW 
0, carbon and oxygen have been ignited and the CO core will 
become an ONeMg core. As described in the paragraph above, 
the star will become a helium giant star after the ejection of its 
hydrogen-rich envelope. If the helium envelope can not get to 
the point of AW = 0, the mass of the ONeMg core might arrive 
at the Chandrasekhar mass limit after about 1000 - 10000 yrs 
and the star might explode as a type lb supernova. Otherwise, an 
ONeMg WD would be its final fate. 

For the cases of extremely low and high metallicities, the sit- 
uation of AW = can not be met at the low-mass end and we 
can not obtain the final mass through AGB stars (see FigQ}. This 
is because for the stars with these extreme condition, the final 
mass is too large with respect to the initial mass and the enve- 
lope is too thin. For example, the mass of the core through AGB 
is larger than 0.6 M for a star of M; = O.8M with Z = 0.0001, 
and the situation of AW = is not obtained before Eggleton's 
stellar evolution code can not work. This may imply that if the 
low-mass stars with extremely low and high metallicities expe- 
rienced the AGB stage and became CO WDs finally, they might 
have a different mechanism to lose their envelopes from those 
with middle metallicities. 



3.2. Hertzsprung-Russel diagram 

The final positions (AW = 0) of the stars before envelope loss 
with the different metallicities in a Hertzsprung-Russel diagram 
(HRD) are shown in Fig. [3] The dotted line gives the division of 
the FGB and AGB stars (the upper populations are AGB stars, 
while the lower ones are FGB stars). The effective temperature 
strongly depends on the metallicity since all the stars are near 
the Hayashi line which is determined by metallicity. 



3.3. comparison between observations and theoretical 
predictions 

A comparison between observations and our theoretical predic- 
tions about IFMRs is shown in Fig. |4] Some empirical rela- 
tions are also presented in the figure. The cross represents the 
mean error of the observational data. The solid line is the IFMR 
for stars of Z = 0.02 obtained in this paper. The dashed, dot- 
dashed and thick dotted lines are the empirical relations given by 
Weidemann (EOOOi l. Ferrario et al. J20051 > and Williams (120061 . 
respectively. It is obvious that our theoretical IFMR for stars of 
Z = 0.02 is consistent with the observations and the empirical 
relations, especially for low mass stars. This is a natural result 
since almost all clusters observed are open clusters in the Galaxy 
and their metallicities are around 0.02. 



1 Note that the mass threshold may have an uncertainty of 0.1 5M Q 
because of the grid density of the models in this paper. The mass interval 
is about 0.2 - 0.3M o for our models. 

2 The non-monotonic effect of metallicity on stellar evolution is from 
equation (1). Please see the fourth paragraph in subsection l3.ll 



3.4. the potential evidence of the dependence of the 
initial-final mass relationship on the metallicities 

The IFMR derived from NGC 2099 might be evidence of the 
dependence of the IFMR on the metallicity. Kalirai et al. (2005 ) 
noticed that half of their data points reside in a region of the 
IFMR s for stars of Z = 0.008 and Z = 0.02 given by Marigo 
(2001 ). Then, they suggested that they would discover the first 
evidence of the effect of the metallicity on the IFMR. We rehan- 
dle the data in Kalirai et al. (2005). The result is shown in Fig[5] 
The solid line and the dot-dashed line are the results in this pa- 
per for Z = 0.02 and Z = 0.01, respectively. The dashed line is 
the best-fit linear least squares line from the data in Kalirai et al. 
(2005). For convenience to compare with the results in this pa- 
per and with the empirical relations, we set the slope of the fitted 
line as that of the empirical relation given by Williams (2006), 
which is consistent with our IFMR for stars of Z = 0.02 (see the 
linear thick dotted line in FigHJ). The method used here is sim- 
ilar to that in calculating a mean value since the mean value is 
the intercept of a line whose slope is zero. When fitting the lin- 
ear line, we use the errors of the data as weight. Seen from Fig 
[5] the fitted line is more consistent with the line of Z = 0.01 than 
that of Z = 0.02. Comparing with the empirical relation given 
by Williams (2006), the fitted line moves upward by 0.05 M . 
Interestingly, O.O5M is equal to the estimate of the AMf be- 
tween the IFMR for stars of Z = 0.02 and that of Z = 0.01 by 
Williams (2006 ). As shown in subsection !3.3l the IFMR for stars 
of Z = 0.02 in this paper is consistent with the observations in 
the Galaxy. Then, although the observational error is large, we 
reconfirm the discovery of Kalirai et al. ((2005) that the IFMR 
derived from observation in NGC 2099 is potential evidence of 
the dependence of the IFMR on the metallicities and that the 
metallicity of NGC 2099 may be about 0.01. Note that the ob- 
servational metallicity of NGC 2099 i s [Fe/ H]~- 0.1, which is 
less than solar metallicity (Kalirai et al. |2005l ). 

Another source of potential evidence is the Hyades. The 
metallicity of the Hyades ([Fe/H]=+0.17) is larger than that 
of NGC 2099, while its age is s imilar to tha t of NGC 2099 
(Perryman et al. |T998l Kalirai et al. 120011120031 ). The mean mass 
of WD in NGC 2099 is < M >= 0.80 + 0.03M o whereas WDs 
in the Hyades have a mean mass of < M >= 0.72 + O.O2M 
(Claver et al. 120011 Kalirai et al. 2005). This also qualitatively 
matches with the trend obtain in this paper that for Z < 0.04 - 
low metallicity, on average, leads to higher final mass. 

Recently, Kalirai et al (12007b! ) showed the first constraints 
on IFMRs at the low mass end. Their results are derived from 
three different clusters with different ages and different metallic- 
ities, i.e. NGC 7789, NGC 6819 and NGC 6791. Since the ages 
and the initial masses for these clusters are different, their results 
can not give a direct constraint on whether the IFMR depends 
on metallicity. However, our results are consistent with that of 
Kalirai et al d2007bl ) within the errors. For example, WD7 in 
NGC 6791 ([Fe/H] = +0.4 + 0.1]) has Mi = U6+°;*}AT and 
M { = 0.53 + O.O2M , which is similar to our results of Z = 0.04, 
i.e. M f = O.5468M for M { = 1.01M o and M f = 0.5560M o for 
M; = 1.31M . The cases of the other two cluster are similar. 
So, Kalirai et al (2007b ) might provide indirect evidence of the 
dependence of the IFMRs on metallicity. 

3.5. the mass distribution of white dwarfs 

Incorporating the IFMR in this paper into Hurley's single stellar 
population synthesis code (Hurley et al. 2000 2002), we simu- 
late the distribution of the final mass. The results for age < 15 
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Fig. 3. The final positions (AW = 0) of FGB/AGB in Hertzsprung-Russel diagram (HRD) for different metallicities. In the upper 
panel, the final positions of AGB stars locate above the dotted line, while the final positions of FGB stars locate under it. The upper 
panel shows models whose metallicity is high enough to permit the low-mass stars to get to the point of AW = on the FGB. Note 
that the scale of the abscissa in the upper panel is different from that in the lower panel. 

Table 1. The masses (in M ) of helium WDs for different metallicities (in Z ) (Row 1) and different initial masses (Column l)(in 
M e ). 



Mi 


1.0 


1.5 


2.0 


2.5 


3.0 


4.0 


5.0 


0.79 
1.00 


0.4554 
0.4697 


0.4415 
0.4518 


0.4311 
0.4463 


0.4257 
0.4354 


0.4185 
0.4318 


0.4108 
0.4213 


0.4074 
0.4176 



Table 2. The final radius (in logRf/R Q ) when AW = in first giant branch for different metallicities (in Z Q ) (Row 1) and different 
initial masses (Column l)(in M Q ). 



Mi 


1.0 


1.5 


2.0 


2.5 


3.0 


4.0 


5.0 


0.79 
1.00 


2.2323 
2.2507 


2.2419 
2.2531 


2.2462 
2.2595 


2.2523 
2.2610 


2.2596 
2.2658 


2.2688 
2.2734 


2.2796 
2.2784 
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Table 3. The maximum masses (Row 3) of CO WDs for different metallicities (Row 1). The initial masses to get the final maximum 
CO WD are shown in Row 2. 



z/z Q 


0.005 


0.015 


0.05 


0.2 


0.5 


1.0 


MJMq 


4.52 


4.32 


4.52 


5.01 


5.50 


5.50 




1.0918 


1.0667 


1.0912 


1.1251 


1.1461 


1.0566 


z/z e 


1.5 


2.0 


2.5 


3.0 


4.0 


5.0 


MJMq 


6.03 


6.03 


6.03 


5.50 


5.50 


4.52 


M™/M Q 


1.1228 


1.1081 


1.1370 


1.0842 


1.2431 


1.2132 



Table 4. Mass threshold (Row 2, in solar mass) for different metallicities (Row 1, in solar metallicity). If a star has an initial mass 
larger than the mass threshold for a certain metallicity, the second dredge-up occurs in the star when helium is exhausted in the 
center. 



z/z 


0.005 


0.015 


0.05 


0.2 


0.5 


1.0 


MJM e 


2.30 


2.30 


2.50 


3.00 


3.30 


3.50 


Z/Zq 


1.5 


2.0 


2.5 


3.0 


4.0 


5.0 


MJMq 


3.50 


3.50 


3.50 


3.30 


3.00 


2.30 




2 4 6 



Fig. 4. Comparison between our theoretical prediction and ob- 
servations. The cross represents the mean error of the obser- 
vational value. Most of the observational data are from the ob- 
servation of some open clusters in the Galaxy and the data are 
from Reid ( 119961 ), Herw ig d 19951 1, Koes ter & Reimers (119961 1, 
Finley & Koester (fl997l>, C laver et al. d200Tb Williams et al. 
( 120041 ), Ferrario et al. (|20D5T > and Dobbie (2006b). The solid line 
is the IFMR for stars of Z = 0.02 in the paper. The dashed, dot- 
dashed and thick dotted lines are the empirical relations given by 
Weidemann (|2000T >, Ferrario et al. ( T2003T ) and Williams (2006), 
respectively. The two thin dotted lines show the range of error of 
the thick dotted line. 

Gyr are shown in Fig. [6] ForZ = 0.02, a constant star formation 
rate of 3M Q /yr is assumed to match the observation in the Galaxy 
(Miller & Scalo [19791 Timmes et al. H997l . For Z = 0.001 and 
Z = 0.0001, a single star burst is assumed. Actually, the star for- 
mation rate does not significantly affect the distribution of the 
final mass. Kepler et al. (2007 ) showed that the mass distribu- 
tion of the single and non-magnetic DA stars in Data Release 4 
of the Sloan Digital Sky Survey (hereafter Sloan DR4) can be 
fitted by four Gaussian fits. The Gaussian fits of the mass dis- 




2 4 6 



Fig. 5. Comparison between our results and the observation in 
NGC 2099 (M37). The cross represents the mean error of ob- 
servational data. The solid line and the dot-dashed one are the 
results in this paper for Z=0.02 and Z=0.01, respectively and the 
dashed one is the best-fit linear least squares line of the observa- 
tional data. The data are from Kalirai et al. (2005 ). 



tribution of DA stars hotter than T eS =12000K from Sloan DR4 
are also shown in Fig. [6] by dotted lines (Kepler et al. 120071 ). 
Incorporating the IFMR for stars of Z = 0.02 into Hurley's rapid 
stellar evolution code, we reproduce the observational mass dis- 
tribution of the WDs except for a small bump at M core ot 1 .2OM . 
The BPS results of Han ( 1998 ) showed that the final mass of the 
merger of two CO WD concentrates in Mf 1 .2OM if the mass 
of the merger is less than the Chandrasekha mass limit (see in 
subsection l4.2l for details). The bump in Sloan DR4 may be de- 
rived from the mergers of double degenerate systems. 

In Sloan DR4, both DA WDs and DB WDs show a trend that 
the mean mass decreases with the effective temperature (Kepler 
et al. 120071 ). Kepler et al. (120071 ) suggested that this trend would 
not be a real case and it would be derived from the absence of 
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Fig. 6. Mass distribution of white dwarfs calculated by Hurley's 
single stellar population synthesis code incorporating our initial- 
final mass relationships for stars of different metallicities. Solid 
line represents the mass distribution of white dwarfs for Z = 
0.02, assuming a constant star formation rate of 3M /yr. Dashed 
and dot-dashed lines represent the mass distribution of white 
dwarfs for Z = 0.001 and Z = 0.0001 with a single star burst, re- 
spectively. The age when WDs form in this diagram is less than 
15 Gyr for all the metallicity. The dotted lines are the Gaussian 
fits of the mass distribution of single and non-magnetic DA stars 
with r eff > 12000K in Dat a Release 4 of the Sloan Digital Sky 
Survey (Kepler et al. 12007] !. 




Fig. 7. Mass distribution of white dwarfs with observational se- 
lection effect. A single starburst is assumed for both metallici- 
ties. The age for every line is shown on the top of the line. 



neutral particles in their model used to fit the WD spectral of 
Sloan DR4. From our study, however, it is very likely that the 
trend is derived from the the effect of the metallicity as suggested 
by Wilson ( 120001 1. We know that stars with a low metallicity 
evolve faster than those with a high metallicity. So WDs from 
stars with a low metallicity form earlier and they have a longer 
time to cool. Meanwhile, WDs with a low metallicity are more 
massive for a certain initial mass if Z < 0.04 (as shown in Figs 
[U and |6]l, and will cool more quickly (Kippenhahn & Weigert 
1990). Both of these facts imply that some massive WDs will 
be cooled and may not be observed. For the reasons above, we 
may expect a higher mean mass of WDs in globular clusters than 
that in open clusters or in the field. However, this prediction was 



not found to hold in NGC 6752, which is an old low-metallicity 
globular cluster (Moehler et al. [2004). The reasons are as fol- 
lows. Most globular clusters are very old and high mass WD 
cools faster than low mass WDs (Kippenhahn & Weigert 1990). 
So, the luminosity of high mass WDs are probably too low to 
be detected. We construct a toy model to examine the validity 
of the explanation. The basic parameters for the toy model are 
from Moehler et al. (2004), i.e. the distance modulus is set to (m- 
M)=13 m .05-13 m .20, which is used to deduce the mean value of 
the WD mass in the globular cluster (0.53- 0.59 M ) by Moehler 
et al. (120041 ); the absolute magnitude of the Sun is 4 m .75; the 
limiting magnitude is set to 26 m since the signal-to-noise ra- 
tio (S/N) is low in Moehler et al. (2004) (see their Fig. 1). The 
lowest luminosity WDs which can be detected in NGC 6752 is 
IQ-i -28 _ io 3 ' 22 L Q , which means that a WD can not be detected 
in the cluster NGC 6752 if it cools for about lGyr (Fontaine 
et al. 120011 Moroni & Straniero 120021 l2007l ). The mass distri- 
butions of WDs with observational selection effect for different 
metallicities and different age are shown in Fig. [7] The initial 
mass function (IMF) is set to <p(m) oc m~ 235 (Salpeter 1955 ) and 
the ages of the progenitors of WDs are estimated by the life of 
the progenitors in main sequence multiplying by a factor of 1 . 1 
(Hurley et al. 2000). The solid and dashed lines show the mass 
distribution of WDs for Z = 0.001 and Z = 0.0001, respectively. 
The age for every line is shown on the top of the line. It is evi- 
dent that the position of the peak in the mass distribution of the 
WDs moves to a lower mass with time. If the age of a globular 
cluster with Z = 0.0001 is larger than about 7Gyr, the luminos- 
ity of all the WDs with M\ > lM will be too low to be detected. 
NGC 6752 is a very old low-metallicity globular cluster (age is 
larger than 12Gyr and z = 0.002, Cassisi et al. 120001 ). For the 
case of Z = 0.001, the peak mass of WDs is smaller than about 
O.6OM at 12Gyr, which is slightly larger than the observational 
value, O.53M - O.59M . Considering an observational error of 
0.05M o and the real metallicity of NGC 6752, our results match 
with the observations of Moehler et al. (2004). A fact should be 
noticed that since NGC 6752 is a globular cluster, the interaction 
among stars in the cluster might reduce the final mass. 

The results of Fig Q indicate that the mass distributions of 
WDs are possibly a good cosmochronometers as first suggested 
by Schmidt (1959 ). On the other hand, the mass distribution of 
WDs also has a potential ability to distinguish the metallicitiy of 
a cluster since the mass peak with a given age is different for dif- 
ferent metallicities. So, the mass distribution of WDs in globular 
clusters can help to break the degeneracy of the metallicity and 
the age if there are enough WDs to be observed. 



4. Discussions and Conclusions 

4.1. the uncertainty of the initial-final mass relationship 

The major uncertainty of the IFMR is from the assumption that 
the final mass is equal to the core mass when the binding en- 
ergy of the envelope AW = for an AGB/FGB star. Although 
the consistency between the mass distribution of WDs in the pa- 
per and that from Sloan DR4 upholds the assumption, there is 
no direct observational evidence to verify it. Another widely ac- 
cepted choice of the final mass is the core mass of a star at the 
beginning of the thermal-pulse AGB (TPAGB, the first TP or the 
end of early AGB) (see the review by Weidemann 2000). From 
this choice, taking Z = 0.02 as an example, the general trend of 
the IFMR is divided into three segments as follows: (i), the core 
mass is almost constant from 1 M to 2.5 M , around 0.52 M , 
which is smaller than that in this paper by about 0.03-0.1 M . 
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Fig. 8. Binding energy of the envelope (panels 1 and 2), helium 
luminosity (panel 3) and stellar radius (panel 4, solid line) as a 
function of time for a star with M\ = 3.00M G and Z = 0.02. 
The time when AW = is set to be zero in the figure. Panels 2 
and 4 are the enlarged images of panels 1 and 3 along abscissa. 
The dependent variable of AW in panels 1 and 2 is the same as 
that in Fig. [2] For comparison, helium and hydrogen luminosi- 
ties are also plotted in panel 4 by dot-dashed and dashed lines, 
respectively, and their values are shown by the ordinate of panel 
3. 



(ii), the core mass strongly increases with the initial mass up to 
4.0 M e . This trend is similar to that in this paper, but the core 
mass is smaller than that in this paper by 0.03-0.05 M . (iii), the 
core mass slowly increases with the initial mass up to 6-7 M . 
The final mass is smaller than that in this paper by 0.1-0.2 M Q . 
According to above discussion , the maximum variation of the 
final mass derived from different assumptions may be as large as 
0.2 M Q . 

In this paper, we skipped the detailed evolution of TP and 
simply treated it as an average evolution, which may lose some 
information on the structural change of the envelope caused by 
thermal pulse at the early phase of thermally pulsing AGB. We 
therefore chose a model with M; = 3.00M o and Z = 0.02 to 
examine the influence of thermal pulses on AW and the result 
is shown in Fig. [8J From the figure, we see that there is a time 
that AW changes its sign but returns back immediately during a 
full amplitude of TP. The change of the sign of AW is from the 
decrease of stellar radius which results from the decrease of to- 
tal luminosity (see panels 3 and 4). This implies that superwind 
might be periodically interrupted, but the timescale of the inter- 
ruption is very short (~ lOOyr), and then its influence is small. 
There are many studies focusing on the TPs previously and we 
briefly summarize them in the following. Generally, the final 
mass after several TPs is larger than the core mass at the first TP. 
Forestini & Charbonnel d 19971 ) calculated a model of Mi = 3M 
assuming a Reimers mass-loss law with 77 increasing from 2.5 to 
5.0 in the final AGB stage. They obtained M f = O.6OM after 
19 TPs, but the core mass is O.54M at the first TP. Dominguez 
( 1999 ) showed that the thermal pulse starts when the core mass 
is 0.57 \Mq for a star of M\ = 3M , while the final mass in- 
creases to Mf = 0.7 1M after 26 TPs if the mass loss prescrip- 
tion of Groenewegen & de Jong ( 1994) was adopted (Straniero 
et al. 119971 ). If the fairly strong mass-loss law of Blocker ( 1995 ) 
was adopted, the final mass becomes O.63M after 20 TPs from 
O.53M at the first TP (Blocker [19951 ) . However, the final mass 



of 0.68 M is obtained from the model with M, = 3.OM and 
Z = 0.02 in this paper. 

In Fig. [8] we see that AW = for a star with Mi = 3.OOM 
and Z = 0.02 is achieved before TPAGB phase, which might be 
very important for galactic chemical evolution since it is widely 
believed that s-process elements are created in TPAGB phase of 
a low-mass star (Han et al. 1995; Busso et al. 119991 ). Actually, 
AW = for all the models in this paper is achieved at EAGB 
phase. These results seem to indicate that the s-process elements 
could not be created. In fact, as mentioned in subsection |2.1| it 
is more likely that AW = is only a lower time limit for super- 
wind, i.e. a superwind starts at or soon after AW = 0. Therefore, 
the stars achieving AW = at EAGB phase still have oppor- 
tunities to enter into TPAGB phase, and contribute to s-process 
elements. So, the final mass shown in this paper might be differ- 
ent from a real one. As mentioned above, the core mass in 3 M 
models grows by about 0.1 M during TPAGB evolution. For 
a 1 .5 M Q star with Z = 0.02, the final core mass is about 0.03 
M & larger than that in this paper if the mass loss prescription 
of Groenewegen & de Jong (1994) was adopted (Dominguez 
1999). The main result in this paper then still holds. 

In Fig. [U we see that when M\ < 3.5M , the IFMR in this 
paper is well consistent with the empirical relations. However, 
for the higher initial mass, the IFMR in this paper is larger than 
the empirical relations given by Weidemann (2000) and Ferrario 
et al. (120051 ). This discrepancy is mainly from our assumption. 
As mentioned in subsection 12. 11 and above paragraph, AW > 
is a necessary condition to eject the envelope of an AGB star. 
Since AW = is achieved at the second dredge-up for the star 
with M\ > 3.5M during EAGB phase and the core is decreas- 
ing at this phase, the necessary condition means that final mass 
may be overestimated. The uncertainty of Mf attributed to the 
definition of core in this paper may also reduce the discrepancy. 
Meanwhile, as mentioned in subsection l3.11 the final mass might 
be overestimated if the helium envelope of a star reaches AW = 
again after the ejection of hydrogen-rich envelope. This overes- 
timate is more likely for high initial mass, since the helium en- 
velope is thicker, and then reaches AW = more easily (see sub- 
sec tion [3~Tl ). However, since the errors from the observations are 
also very large (see the error bar in Fig. [4|, the IFMR obtained 
in this paper is well located in the error range and therefore, is 
still consistent with observations. 



4.2. low-mass white dwarf 



Kalirai et al. d2007al ) suggested that helium white dwarfs may 
be obtained in metal-rich old clusters, i.e. NGC 6791. It is a 
reasonable assumption that the same mechanism in NGC 6791 
would work for metal-rich field stars. Recently, Kilic, Stanek & 
Pinsonneault d2007d l rehandled the data in Valenti & Fischer 
(2005) and found that there have been metal -rich stars with 
[Fe/H]> at all times in the local Galactic disk, although the 
metallicity distribution of disk stars peaks below solar metallic- 
ity for stars with ages greater than about 5 Gyrs. Considering that 
only 5% of all WDs in the local disk are single low-mass white 
dwarfs (< 0.45M o ) (Liebert et al. 120531 Kepler et al. 12007] ) and 
the fraction of metal-rich stars with ages greater than 9 Gyrs is 
21%, Kilic, Stanek & Pinsonneault (2007c) argued that only the 
stars of [Fe/H]>+0.3 can lose their hydrogen-rich envelope to 
produce helium WDs on the FGB. From our study, however, all 
stars with M\ < 1 .OM and Z > 0.02 will lose their hydrogen- 
rich envelope and finally become helium white dwarfs. This in- 
consistency may result from the overestimate of metal-rich stars 
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in the sample of Valenti-Fischer (Reid et al. 120071 1 adopted by 
Kilic, Stanek & Pinsonneault (|2007cl > 

As shown in our study, there should be numerous He WDs in 
metal-rich old clusters, which has indeed been observed in NGC 
6791 (Kalirai et al. 2007a). Here, we emphasize that we do not 
assume any mass-loss mechanism while Kalirai et al. (2007a) 
assumed a metal-enhanced stellar wind on the red giant branch 
(RGB). Our study shows that AW = in several FGB stars with 
high metallicities and small initial masses can be achieved before 
these stars get to the tip of FGB, indicting that the relative num- 
ber of the tip RGB stars in metal-rich old clusters is smaller than 
that in metal-poor clusters. This fact is directly observed in the 
RGB luminosity functions of two open clusters, e.g. NGC 188 
and NGC 6791 (see Fig 8 in Kalirai et al |2007al l. Meanwhile, 
a similar effect should be seen in the local population of RGB 
stars. Luck & Heiter (see Fig 9 in their paper, 2007) make a 
comparison of the metallicity histograms between dwarfs and 
giants within 15 pc of the Sun and found that the dwarf pop- 
ulation has a high metallicity tail extending up to [Fe/H]~ 0.6, 
while the giants show a sharp drop in numbers after [Fe/H]=0.2 
and no giants with [Fe/H]> 0.45 are observed in the field. 

The fate that low-mass metal-rich stars will become He WDs 
before helium is ignited on FGB might give some constrain on 
planetary nebulae (PNe). It is widely accepted that PNe originate 
from AGB stars or are related to binary evolution. For the case 
of binary evolution, one component in the binary system fills 
its Roche lobe at AGB phase and the mass transfer is dynami- 
cally unstable which leads to a common envelope phase. After 
the ejection of the common envelope, a PN may form. However, 
since some metal-rich stars may not experience AGB phase, we 
might speculate that the number of PNe in metal-rich galaxies 
is relatively lower than that in metal-poor galaxies. Interestingly, 
Buzzoni et al (2006, see their Fig 11) really observed a rela- 
tively low number of PNe per unit galaxy luminosity in more 
metal-rich elliptical galaxies. Gesicki & Zijlstra ( 120071 ) com- 
pared the mass distribution of the central stars of planetary neb- 
ulae (CSPN) with those of WDs. These two distribution are very 
different, i.e. the CSPN mass distribution is sharply peaked at 
0.61 M Q ranging from 0.55 M Q to 0.66 M Q , while the WD distri- 
bution peaks at a slightly lower mass and shows a much broader 
range of masses. Gesicki & Zijlstra (2007) suggested that this 
difference may imply that only some WDs have gone through 
the PN phase. Our models provide a channel for WDs to avoid 
the PN phase. 

Based on our results, WDs from single stars are always larger 
than 0.4 M Q , whatever the composition of the WD is. Then, 
extremely low-mass WDs ( O.2M ) are possible in binary sys- 
tems or once in binary systems since no stellar population is 
old enough to produce such extremely low-mass WDs through 
single star evolution. Only recently, several extremely low-mass 
( 0.2M o ) W Ds were discovered in th e field (Liebert et al. 120041 
Kawka et al. 120061 Eisenstein et al. 120061 Kilic et al. I2007al) . 
and some of them are companions of pulsars (van Kerkwijk et 
al. 120051 Kilic et al. d2007bl ) also found a low-mass WD in a bi- 
nary system. If an extremely low-mass WD were a single star, it 
is very likely that it could have a very high space velocity since 
it may come from a close double-degenerate binary, where its 
companion has gone through a supernova event that disrupted 
the binary (Hansen 2003). LP 400-22 might be a case from this 
channel (Kawka et al. 12006) . In any case, hard work is need to 
systematically search for the companions of WDs with mass less 
than 0.4 M . 



4.3. the potential effect of the initial-final mass relationship 
on type la supernovae 

As the best cosmological distance indicators, Type la supernovae 
(SNe la) have been successfully applied to determine the cosmo- 
logical parameters ,e.g. Om and £2a (Riess et al.[T998; Perlmutter 
et al. 1999). Phillips relation (Philips 1993 ) is used when taking 
SNe la as the distance indicators. It is assumed that Phillips re- 
lation is correct at high redshift, although the relation was ob- 
tained from a low-redshift sample. This assumption is precari- 
ous since the exact nature of SNe la is still unclear, especially 
the progenitor model and explosion mechanism (Hillebrandt & 
Niemever 120001 Leibundgut 2000). If the properties of SNe la 
evolve with redshift, the results for cosmology might be dif- 
ferent. Since metallicity decreases with redshift, a good way to 
study the correlation between the properties of SN la and redshift 
is to study the correlation between the properties of SN la and 
metallicity. It is widely believed that SNe la are from thermonu- 
clear runaway of a carbon-oxygen white dwarf (CO WD) in a bi- 
nary system. The CO WD accretes material from its companion 
to increase its mass. When its mass reaches its maximum stable 
mass, it explodes as a thermonuclear runaway and almost half 
of the WD mass is converted into radioactive nickel-56 (Branch 
2004). The mass of nickel-56 determines the maximum lumi- 
nosity of SNe la. The higher the mass of nickel-56 is, the higher 
the maximum luminosity is (Arnett 1982). Some numerical and 
synthetical results showed that metallcity may affect the final 
amount of nickel-56, and thus the maximum luminosity of SNe 
la (Timmes et al. 120031 Travaglio et al. 120051 Podsiadlowski et 
al. 12006 ). There is also much evidence about the correlation be- 
tween the properties of SNe la and metallicity in observations 
(Branch & Berg h [19931 Hamuy et al. [19961 Wang et al. [19971 
Cappellaro et al. 1 19971 ). e.g. the maximum luminosity of SNE la 
is proportional to the metallicity (Shanks et al. 2002 ). 

Two progenitor models of SNe la have competed for about 
three decades. One is a single degenerate model, which is 
widely accepted (Whelan & Iben 119731 . In this model, a CO 
WD increases its mass by accreting hydrogen- or helium-rich 
matter from its companion, and explodes when its mass ap- 
proaches the Chandrasekhar mass limit. The companion may be 
a main-sequence star (WD+MS) or a red-giant star (WD+RG) 
(Yunge lson et al. [19931 Li e t al. [19971 Ha chisu et al. I1999al 
[1999b j Nomoto et a l. [19991 Langer et al. 120001 . Hachisu & 
Kato d2003al I2003bt suggested that supersoft X-ray sources, 
which belong to WD+MS channel, may be good candidates 
for the progenitors of SNe la. The discovery of the companion 
of Tycho's supernova also verified the reliability of the model 
(Ruiz-Lapuente et al. |2004t Ihara et al. 2007 ). Another progen- 
itor model of SNe la is a double degenerate model (Iben & 
Tutukov 1984; Webbink 119841 ), in which a system consisting 
of two CO WDs loses orbital angular momentum by gravita- 
tional wave radiation and merges. The merger may explode if 
the total mass of the system exceeds the Chandrasekhar mass 
limit (see the reviews by Hillebrandt & Niemeyer 2000] and 
Leibundgut 120001 ). In both of the progenitor models, the CO 
WD which finally explodes as a SN la should approach or ex- 
ceed the Chandrasekhar mass limit. Obviously, a higher mass 
CO WD may fulfill this situation more easily than a low mass 
one. According to our results, the mass of a CO WD with a 
given initial mass will be higher in the circumstance with ex- 
tremely high metallicity or extremely low metallicity than that 
in the middle-metallicity circumstance. Metallicity is therefore 
very relevant to SNe la. 
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4.4. Conclusions 

We use the method of Han et al. (119941 1 to calculate the IFMR for 
stars of different metallicities. The conclusions are as follows. 

1 . There is an obvious dependence of the IFMRs on the metal- 
licity and the dependence is not monotone. When Z = 0.04, 
the final mass of the CO WD for a given initial mass is small- 
est. For higher or lower Z, the mass of CO WD will be higher. 
The difference of the final mass derived from different metal- 
licties is up to O.4M . 

2. The initial-final mass relationship for stars of Z = 0.02 is 
consistent with observations. 

3. For Z > 0.02, a helium white dwarf is formed from a star of 
Mi < I.OMq. The final masses of helium WDs for a given 
initial mass slightly decreases with metallicity. 

4. Incorporating the IFMR for stars of Z = 0.02 in the paper 
into Hurley's single stellar population synthesis code, we re- 
produce the mass distribution of DA WDs in Sloan DR4 ex- 
cept for some extra-massive WDs. 

5. We reconfirm the discovery of Kalirai et al. (2005) that 
the initial-final mass relationship derived from observation 
in NGC 2099 might be evidence of the dependence of the 
IFMR on metallicities and that the metallicity of NGC 2099 
may be about 0.01, although the observational error of white 
dwarfs in NGC 2099 is large. It should be encouraged to pro- 
gram more accurate observations to find a larger sample of 
white dwarfs in globular clusters. Such programs may help 
to confirm the dependence of the initial-final mass relation- 
ship on the metallicity. 

6. We bring up again Willson's suggestion (Willson 2000) that 
the effect of metallicity may be the origin of the phenomenon 
that the mean mass of WDs decrease with effective tempera- 
ture. 
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Appendix A: FITTED FORMULAE OF INITIAL-FINAL 
MASS RELATIONSHIP FOR STARS of 
DIFFERENT METALLICITIES 

The initial-final mass relationship for stars of different metallic- 
ities can be approximated by two parabola segments except for 
the case of Z = 0.0001. In this Appendix, we present the fitted 
formulae. Hereafter, M\ is the initial mass of stars and Mf is the 
final mass of stars. 
ForZ = 0.1, 

M f = max[M f ,i,Mf, 2 ], l.OM < M < 5.5M ,err < 1.0%, (A.l) 
where 

M fjl = 0.5675 + 0.1027M; - 0.0006261M 2 , (A.2) 
and 

Mf, 2 = 0.7345 - 0.07362Mi + 0.03973M 2 . (A.3) 
For Z = 0.08, 

M f = minLMf^Mf^], 1.0M o < M { < 6.OM ,err < 1.8%, (A.4) 
where 

Mf.j = 0.7097 - 0.1930Mi + 0.08235M 2 , (A.5) 
and 

Mf, 2 = 0.8316 -0.1 183Mj+ 0.0351 1M 2 . (A.6) 
For Z = 0.06, 

M f = min[M f ,i,Mf, 2 ],O.8M < M ; < 6.5M ,err < 2.3%, (A.7) 
where 

M f ,! = 0.5806 - 0.06852Mi + 0.03928M 2 , (A.8) 
and 

Mf, 2 = 0.8957 - 0.1313Mi + 0.03004M 2 . (A.9) 
For Z = 0.05, 

M f = min[M f ,i,Mf, 2 ],O.8M < M { < 6.5M ,err < 1.8%, (A.10) 
where 

M f ,i = 0.5736 - 0.06446Mi + 0.03621M 2 , (A. 11) 

and 

Mf, 2 = 0.8971 - 0.1255Mi + 0.02750M 2 . (A. 12) 

For Z = 0.04, 

M f = mm[Mf A ,M fa ],0.^M o < M { < 6.5M ,err < 1.9%,(A.13) 
where 



M f ,! = 0.5737 - 0.06207Mi + 0.03353M 2 , (A. 14) 
and 

Mf, 2 = 0.8691 -0.1 107Mi + 0.02491M 2 . (A.15) 
For Z = 0.03, 



M f = min[M f ,i,Mf, 2 ],O.8M < M ; < 6.5M ,err < 2.3%, (A.16) 
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where 

M f ,i = 0.5701 - 0.05225M; + 0.030 13M 2 , (A. 17) 
and 

M f>2 = 0.9897 - 0.1608Mi + 0.03022M 2 . (A. 18) 
For Z = 0.02, 

Af f = min[Af f ,i,M f , 2 ],O.8Af < M { < 6.0M G ,err < 2.1% (A. 19) 
where 

M f ,i = 0.5716 - 0.04633M; + 0.02878M 2 , (A.20) 
and 

M f>2 = 1.1533 - 0.2422Mi + 0.04091M 2 . (A.21) 
ForZ = 0.01, 

Mf = min[M f ,i,M fj2 ],0.8M o < M { < 5.5M ,err < 2.2% (A.22) 
where 

M f>1 = 0.5897 - 0.0563 1M; + 0.03395M 2 , (A.23) 
and 

M f>2 = 0.8660 - 0.1240Mi + 0.03183M 2 . (A.24) 
For Z = 0.004, 

Mf = min[M u ,M f)2 ],0.8M o < M { < 5.0M G ,err < 2.2%, (A.25) 
where 

M tl = 0.5983 - 0.04908M; + 0.03780M 2 , (A.26) 
and 

M f>2 = 0.8529 - 0.1056Mi + 0.03194M 2 . (A.27) 
ForZ = 0.001, 

M f = min[Mf,i,M f , 2 ],0.8M o < M; < 4.5M Q ,err < 1.0%, (A.28) 
where 

M f ,, = 0.5959 - 0.007419Mi + 0.03297M 2 , (A.29) 
and 

M f>2 = 0.6312 - 0.04458M; + 0.04213M 2 . (A.30) 

For Z = 0.0003, 

M f = min[Mf,i,M f , 2 ], 1.0M o < M ; < 4.5M ,err < 0.4%, (A.31) 
where 

M f ,, = 0.6300 + 0.007371Mi + 0.02700M 2 , (A.32) 
and 

M f>2 = 0.7629 - 0.03288Mi + 0.02390M 2 . (A.33) 

ForZ = 0.0001, 

M f = 0.6033 + 0.06839M; + 0.008401M 2 , 

1.0M o < M; <4.5M Q ,err < 1.0%. ( ' 



